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Abstract 

We investigate the Landshoff-Nachtmann two-gluon-exchange model 
of the Pomeron using gluon propagators computed in the Landau gauge 
within quenched lattice QCD calculations. We first determine an effective 
gluon-quark coupling by constraining the Pomeron-quark coupling to its 
phenomenological value (3q = 2GeV _1 . We then provide predictions for a 
variety of diffr active processes. As the propagators have been evaluated 
entirely from QCD first principles (although in the quenched approxima- 
tion), our results provide a consistency check of the Landshoff-Nachtmann 
model. We address the issue of the possible gauge-dependence of our re- 
sults, which will be the object of a future study. 



1 Present address: D.A.M.T.P., University of Liverpool, Liverpool L69 3BX, U.K. 



1 Introduction 



The description of diffractive hadronic physics in terms of the exchange of a 
"Pomeron" has proved remarkably durable, and has successfully withstood the 
advent of QCD. Interest in the Pomeron has been renewed by recent results from 



Total cross sections and diffractive processes are essentially soft in nature, 
and thus do not fall within the realm of perturbative QCD. Thus any fixed-order 
perturbative calculation must be regarded purely as a model of the interaction, 
but one which may be representative of an all-orders QCD result. Nevertheless, 
a successful phenomenological model of Pomeron exchange is provided by two- 
gluon-exchange (2GE) [I], |2j. Crucial to this model is that the gluons be non- 
perturbative, i.e. infra-red finite, both to avoid the Coulomb singularity at t = 0, 
and to provide the correct t-dependence of the differential cross section || ||. 

Recently there has been a series of papers || |6j investigating the phenomenol- 
ogy of the 2GE model in the form proposed by Landshoff and Nachtmann (LN) ||. 
In these papers a non-perturbative gluon propagator extracted from the solution 
of the Schwinger- Dyson equation J7| is inserted in the LN model. Here we will 
explore the phenomenology of the same model, using gluon propagators extracted 
from lattice calculations in the Landau gauge. 

Our method has the advantage that one is inserting in the model a genuine, 
nonperturbative QCD quantity, which has been computed from first principles 
rather than derived from an approximate equation. Also, as we will explain in 
the following, we have only one free parameter in our approach, whereas the 
Schwinger-Dyson solution depends on two parameters. For these reasons, we 
think that our approach provides in principle the best chance to test the success 
of the LN model from the point of view of QCD. 

The layout of the remainder of this paper is as follows. In the next section 
we will summarise the lattice calculations, and present the results for the lattice 
gluon propagator. In section 3 we will describe how we apply the measured gluon 
propagator to the calculation of cross-sections, and compare our results to the 
experimental data. In section 4 we comment on the issue of gauge dependence of 
our results. Finally we present our conclusions. 

2 Computational details 

We use the standard lattice definition of the gluon fields in terms of the link 
variables f|, 
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P 


Size 


No. Cfgs. 


a- 1 (GeV) 


This work 


6.0 


16 4 


150 


1.9 


This work 


6.2 


16 4 


150 


2.7 


Reference lfj 


6.0 


24 3 x 48 


500 


1.9 



Table 1: Parameters of the lattices used in the calculation. The quoted value of 
the lattice spacing is obtained from an analysis of the string tension. 



where a is the lattice spacing and go the bare coupling constant. The gauge con- 
figurations are fixed to the Landau gauge by imposing the gauge-fixing condition 

A(x) = £ A^x + e„) - A^x) = 0. (2) 

The accurate implementation of this step is crucial, and we use the Fourier- 
accelerated algorithm of reference 0. 

The unrenormalised gluon propagator in momentum space is obtained from 
the gauge-fixed gluon fields by taking the Fourier transform 

D, u (p) = (AMM-P))- ( 3 ) 
In the continuum limit, we can write the Landau-gauge propagator in the form 

^(p) = |^-^}Aat(p). (4) 

We performed calculations on two sets of hypercubic lattices, corresponding 
to different physical volumes and lattice spacings. In addition, we used data for 
the gluon propagator as evaluated in Reference [|IJ; such data are particularly 
useful for an analysis of finite volume effects in our calculations. The parameters 
used are listed in Table |l|. Full details of the method and of the simulation for 
the hypercubic lattices are contained in ref. | IT| . 

As is customary in lattice QCD, we fix the value of the lattice spacing in units 
of energy from string tension measurements. This sets the scale of momenta for 
the gluon propagator. 

A representative propagator, D iat (p), is shown in Figure ([!]). It is important 
to keep in mind that by defining the theory on a finite lattice we introduce both 
an ultraviolet and an infrared cutoff for Greens functions. For this reason, our 
gluon propagator is infrared finite by construction. On the other hand, when 
the lattice momenta are such that \pa\ > 1, the theory becomes sensitive to the 
ultraviolet cutoff and one has to worry about finite lattice spacing effects. It is 
therefore crucial to probe the sensitivity of our results to the values of the cutoff 
scales by changing the lattice parameters as discussed above. 
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Figure 1: The gluon propagator on a 16 4 lattice at (5 = 6.0. 
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Size 


a 


a\ 


a 2 


b 


h 


6.0 


16 4 


3.37(1) 


0.176(1) 


0.445(5) 


3.139(4) 


0.193(1) 


6.2 


16 4 


2.72(1) 


0.122(1) 


0.461(5) 


2.645(4) 


0.133(1) 



Table 2: Parameters of the fit to the gluon propagator on the hypercubic lattices, 
in units of the lattice spacing, as discussed in the text 



We fit the propagator to the form 

Aat(p) = 



bo/(p 2l+b2 ) 



\pa\ < 1 
\pa\ > 1 



(5) 



The parameters of the fits in lattice units for each of our hypercubic lattices are 
shown in Table for the remaining lattice, we use the parameters quoted in 
Reference [ffUl. 



3 Results 
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3.1 Determination of an effective quark-gluon coupling 

The salient feature of the LN model is that Pomeron exchange between quarks 
behaves like a C = +1 photon-exchange diagram, with amplitude 

i^(u llx u)(u^u). (6) 

/3q represents the strength of the Pomeron coupling to quarks, and is related to 
the (non-perturbative) gluon propagator by 



P 2 o = ^- 2 jd 2 P ^D( P )\-, (7) 
where g is the gluon-quark coupling. /3 can be determined from, for example, 



2 



the total pp cross section 12], [13 



As a consequence, the LN model yields simple formulae for pp scattering, ex- 
clusive p production in deep inelastic scattering and the J/^— nucleon total cross 
section, which all contain integrals in momentum space of g 2 D(p) 0, In order 
to extract predictions from the model, one needs an expression for g 2 D(p) which 
is valid both in the low momentum region p sa 1 GeV, where in fact the dom- 
inant contribution to the pomeron is expected, and for higher momenta, where 
the perturbative behaviour is recovered. Obviously, convergence of Equation (0) 
requires that the infra-red pole of the gluon propagator be removed by some 
nonperturbative mechanism. 

In order to use the lattice gluon propagator Aat(p) in the LN model we make 
the following assumptions: 

1. we neglect the running of the QCD coupling, i.e. we make the approxima- 
tion g(p) = g; 

2. we assume that in the continuum limit -Di at (p) is multiplicatively renormal- 
isable, as it is in perturbation theory. 

As the scale for the momenta in Di at (p) is set from string tension measurements, 
our assumptions imply that we only have one free parameter to fix in the expres- 
sion g 2 D\ at (p). This is a multiplicative factor, that corresponds to the product 
of a gluon wavefunction renormalisation constant times a numerical value for g 2 . 
We will call this parameter g 2 s , although strictly speaking it is more than just a 
squared coupling constant. 

g 2 s can be determined by using Equation (0) as a normalisation condition, 
i.e. by imposing that j3 attains its phenomenological value of 2.0 GeV -1 : 
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^ = 36^ J d p \- 9eS Aat(p) J =4GeV ( 8 ) 
It is the combination g 2 s D\ at (p) that we insert in the formulae of the LN model. 
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Figure 2: The integral /(s 1 / 2 ) against s 1 / 2 on the 16 4 lattice at (3 = 6.0 

As we mentioned in the previous section, one has to remember that the mo- 
mentum dependence of the lattice gluon propagator is influenced by lattice arti- 
facts. At high momentum, when p becomes comparable to the ultraviolet cutoff 
a" 1 , the propagator is affected by discretisation errors. In order to judge their 
likely importance in this calculation, we show in Figure (0) the integral 



1 



,1/2 



^ 1/2 ) = ^y o <pp[rh*w> (9) 

that we use to constrain g e g, for the 16 4 lattice at (3 — 6.0. The dominant 
contribution to the integral indeed arises from momenta with \pa\ < 1, suggesting 
that our evaluation of the integral Equation (0) should suffer only weakly from 
discretisation uncertainties. A similar conclusion can be drawn from the other 
lattices in this study. 

For low momenta, finite volume effects are present, making it difficult to 
identify a possible gluon mass. In order to isolate the extent to which such 
effects contaminate our calculation, we need a comparison with a calculation at 
a different lattice volume, but at the same lattice spacing. Such a comparison 



is provided by Reference [10], where a high-statistics measurement of the gluon 
propagator is made on a 24 3 x 32 lattice at /3 — 6.0. In the remainder of this 
paper, we will quote results using this determination of the lattice propagator, 
as well as our own determination. 
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Size 


ti (GeV*) 


This work 


6.0 


16 4 


0.85 


This work 


6.2 


16 4 


1.03 


Reference [jlOjj 


6.0 


24 3 x 32 


0.93 



Table 3: The value of /ig defined by Equation (|T0|) on each lattice. 



The inverse "radius" of the Pomeron, /xq, can be determined from the study 
of exclusive p-production in Deep Inelastic Scattering. It is related to the non- 
perturbative gluon propagator by 



2 _S<Ppp 2 [g 2 D(p)Y 



and thus is independent of the coupling g in our approximation, and furthermore 
should not depend on how we choose to parameterise the measured lattice gluon 
propagator. The phenomenological value is liq = 1.1 GeV. The determination of 
Ho on each lattice is shown in Table (P); the consistency between the different 
measurements is encouraging, and the results close to the phenomenological value. 

We now have the ingredients to explore the phenomenology of the 2GE model: 
the p 2 dependence of the propagator determined from a lattice calculation, and 
an effective gluon-quark coupling determined from the fit to (3q. We begin the 
exploration with the analysis of elastic proton-proton scattering. 

3.2 Proton-proton elastic scattering 

The total and elastic differential cross section for proton-proton scattering pro- 
vides a benchmark for the 2GE model of the Pomeron. Indeed it has already 
been explored extensively using nonperturbative propagators obtained from the 
approximate solution of the Schwinger-Dyson equation [|l^, ||. We review the 



formalism here, and then apply the model using the lattice gluon propagator. 



The 2GE amplitude for elastic proton-proton scattering can be written [16], |17 

A(s, t) = is8a 2 s [T x - T 2 ) (11) 



with 



Ti = f d 2 kD(q/2 + k)D(q/2 - k) [G p (q, 0)] 2 (12) 
Jo ' 

T 2 = J" d 2 kD(q/2 + k)D(q/2-k)G p (q,k-q/2) 

[2G p (q,0)-G p (q,k-q/2)} (13) 
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P 


Size 


ff t°ot ( mb ) 


B (GeV^) 


This work 


6.0 


16 4 


18.12 


13.6 


This work 


6.2 


16 4 


19.21 


12.9 


Reference 


6.0 


24 3 x 32 


19.85 


12.6 



Table 4: The energy-independent total cross section, af ot , and the logarithmic 
slope parameter, B. 



where G p (q,k) represents a convolution of proton wave functions. Here T\ and 
T 2 represent the contributions when both gluons attach to the same quark and 
to different quarks within the proton respectively; our expectation [|J, motivated 
by the additive quark rule, is that T\ should dominate T 2 . Following ref. [|15|], we 
take G p (q, k — q/2) = F\(q 2 + \k 2 — g 2 /4| ) , where F\ is the elastic isoscalar form 
factor of the proton, 

4m 2 -2.79t 

l(t) ~ (4m2-t)(l-t/0.71) 2 ' ( } 

The total and elastic differential cross-sections are given by 

< = ^ (15) 

IS 

— = (16) 

dt 167TS 2 

respectively. Single Pomeron exchange is expected to dominate the differential 
cross-section out to —t — 0.5 GeV 2 . To fully describe the energy dependence, 
the intercept of the Pomeron trajectory is taken to be somewhat larger than 1, 
and the measured total and differential cross sections are related to the energy- 
independent expressions above by 

'tat = (i) <?L (17) 

da 
~dl 

For small t, the elastic differential cross section behaves like e Bt , and the model 
is characterised by two parameters, cr° ot and B. 

Following the discussion in the previous subsection, we compute cr° ot and B 
on each lattice using the lattice gluon propagator and the effective coupling g e g. 
Our results are presented in Table |j. The general agreement between the values 
of of ot and B obtained on the various lattices is striking, suggesting either that 
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Figure 3: Data for the pp elastic cross section at ^fs = 53 GeV from Ref. [[jj 
together with the lattice prediction corrected for energy dependence on the 16 4 
lattices at /3 — 6.2 (solid) and (3 = 6.0 (dots), and on the 24 3 x 32 lattice at 
(3 = 6.0 (dashes). 



both quantities are subject to only small discretisation or finite volume effects, 
or that such effects are absorbed in g e g. They are also encouragingly close to the 
phenomenological values of of ot ~ 22.7 mb and B ~ 11 GeV -2 . 

The quality of the fit to the data can be seen in Figure (|3|), where we show 
ISR data for the differential elastic cross section at yfs = 53 GeV together with 
the lattice predictions, with the energy correction of Equation flT8j). Though 
the agreement of the lattice computation with the experimental data is good, 
differences between the lattice results are now clearer, and suggest the need to 
repeat the calculation at different volumes and at different lattice spacings. In 
particular, our lattice at /3 — 6.2 corresponds to a somewhat tiny physical volume 
for hadronic physics, so that finite volume effects may be non-negligible in that 
case. 

3.3 J/ ^-Nucleon Scattering 

This process provides a further important test of the 2GE model, and in the 
following we adopt the analysis procedure of Ducati et al. |]]. The amplitude for 
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Figure 4: The meson-nucleon total cross section as a function of the radius used 
in the meson form factor, Equation (p0[), on the 16 4 lattices at (3 — 6.2 (solid) 
and (3 = 6.0 (dots), and on the 24 3 x 32 lattice at (3 = 6.0 (dashes). Also shown 
are the radii and total cross sections corresponding to the ir, K and J/i/j. 



meson-nucleon scattering is 



A(s,t) = is^a 2 s 



d 2 kD(k)D(2q-k) 



fM{q 2 )-f M {{q-kf) 



f N {q 2 )-fN{q 2 -\q-k + \k 2 ) 



We use the pole approximation for the form factors 

1 



fi{k 2 



(19) 



(20) 



taking (r*> = 0.67 fin 2 , (r£) = 0.44 W, (r 2 K ) = 0.35 W and (rj /ip ) = 
0.04 fm 2 . 

Two phenomenological features are evident in the calculation. Firstly, the 
quark- counting rule is closely satisfied in the case of the hadrons composed of 
light quarks: 

/5 = 6.0, 16 4 : 0^/(7^ = 1.8, 

(3 = 6.2, 16 4 : (Tpp/cTvp = 1-6, 

(3 = 6.0, 24 3 x 32: a pp /a np = 1.5. 



(21) 
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Secondly, the Pomeron couples more weakly to mesons composed of heavier 
quarks. This can be seen in Figure (f|) where we show the meson-nucleon cross 
section as a function of the pole radius used in Equation (pCf), together with 
a Regge fit to the energy-independent part of the n~p and K~p cross sections 
performed in Reference pt| . The differences between the results on the different 
lattices reinforces the need to repeat the calculation on a wider range of lattice 
parameters. 

4 Gauge Dependence 

The issue of gauge dependence in the LN model is a difficult one. Indeed, the 
starting point of LN is the leading-order perturbative calculation of a quark-quark 
scattering amplitude, with the constraint that no colour is exchanged between 
quarks. This is of course a gauge invariant quantity, expressed by a two-gluon 
exchange diagram. On the other hand, LN argue that the salient features of 
pomeron physics can be captured by substituting in such a diagram the tree-level 
gluon propagator with an effective one, including non-perturbative self-energy 
effects, which should in particular account for the removal of the infrared pole. 
Clearly any ansatz for the self-energy introduces in principle an uncontrollable 
gauge dependence in the calculation. Given the phenomenological success of the 
model, one is led to think that maybe the gauge dependence intrinsic in the 
propagator is either mild by itself or gets suppressed (cancels) in the evaluation 
of physical quantities. We emphasise that this problem is inherent in the model 
and has nothing to do with using a lattice propagator. From the point of view 
of our method, one may speculate that the simplest conceivable mechanism for 
the suppression of gauge dependence would be if it gets "factored out" in the 
definition of g e s. In other words, a gauge dependent definition of g e s may cancel 
most of the gauge dependence in the momentum space integrals of the propagator. 

We plan to investigate directly the issue of gauge dependence in a future work 
by using lattice propagators obtained in different gauges. 

5 Conclusions 

We have demonstrated how the Landau-gauge lattice gluon propagator employed 
in the two-gluon-exchange model of proton-proton scattering provides a highly 
successful description of the data, where the only parameter that needs to be 
fixed is an effective quark-gluon coupling. For the first time a genuine QCD 
quantity, evaluated entirely from first principles, has been inserted in the LN 
model, providing an important consistency check of the model itself. Also, given 
the fact that in our approximation the effect of quark loops diagrams is absent, 
our analysis shows the robustness of the LN model, despite its simplicity. 
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Although the description of the data appears largely independent both of the 
volume of the lattice and of the lattice spacing, a more detailed analysis of lattice 
systematic errors is necessary. We plan to perform such an analysis, together with 
an investigation of the gauge dependence of our results, in a future publication. 
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